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Abstract 

We consider a finite or countable collection of one-dimensional Brownian particles whose 
dynamics at any point in time is determined by their rank in the entire particle system. Using 
Transportation Cost Inequalities for stochastic processes we provide uniform fluctuation bounds 
for the ordered particles, their local time of collisions, and various associated statistics over 
intervals of time. For example, such processes, when exponentiated and rescaled, exhibit power 
law decay under stationarity; we derive concentration bounds for the empirical estimates of the 
index of the power law over large intervals of time. A key ingredient in our proofs is a novel upper 
bound on the Lipschitz constant of the Skorokhod map that transforms a multidimensional 
Brownian path to a path which is constrained not to leave the positive orthant. 



1 Introduction 

Define the set I as {1, ... , K} for some K E N or as the set of natural numbers N and let Si, i € I 
be a finite or countable collection of real constants. Consider the following system of stochastic 
differential equations: 

dX l {t) = ^r^{x l (t)=x (]){ t)}dt + dW l {t), iel. (LI) 

Here, Xn\(t) < Xo2)(t) < • • • are the coordinates of the process Xi(t), i G I in the increasing order 
and W = (Wi : i E I) is a system of jointly independent one-dimensional standard Brownian 



motions. The equations in (I.l) model the movement of the particles by interacting Brownian 
motions such that at every time point, if we order the positions of the particles, then the i-ih 
ranked particle from the bottom gets a drift Si for every i E I. As time evolves, the Brownian 
motions switch ranks and drifts, and, hence, their motion is determined by these time-dependent 
interactions. When I is finite, the existence and uniqueness in law of such processes is a consequence 



of Girsanov's Theorem (see e.g. Lemma 6 in 1371). The countable case is subtle and requires 



constraints on the initial positions of the particles. In particular, for (I.l) to make sense, the 
number of particles on every interval of the form (— oo, x] has to be finite at any point in time with 
probability one. We discuss this issue in more detail later in the text. 
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Different versions of the particle system in (1.1) have been considered in several recent articles. 



Among the more recent ones, see Banner, Fernholz and Karatzas [4j, Banner and Ghomrasni [5], 
McKean and Shepp |33|, Pal and Pitman [37J, Jourdain and Malrieu |28j . Chatterjee and Pal 0|8], 
Ichiba and Karatzas [26J, Ichiba et al. [27J, and Shkolnikov [42, 43] . We refer the reader to the above 
articles for the full list of applications of such processes. Related discrete time processes are studied 
in the context of the Sherrington-Kirkpatrick model of spin glasses by Arguin and Aizenman [2], 
Ruzmaikina and Aizenman [JU] and Shkolnikov |41j . 

Classically, the case when <5j = for all i & I has been treated by Harris [23J, Arratia [3] and 
Sznitman |45|.|46| p. 187], and more recently by Swanson [44] . In this case one imagines countably 
many Brownian particles moving under elastic collision, i.e. bouncing off one another when their 



paths meet. This is the process of ordered particles derived from the system in (1.1). Harris |23[ 
(7.1)] gave an explicit formula for the law of Bt \(t), the location at time t of the particle which is 
the leftmost in the infinite system with 5{ = 5 for all i E I. From this he deduced for 25 = 1 that 

-^(O)^) d /2\ 1/4 1 

■ ' 1 : -B(l) as t -)• oo (1.2) 



i 1 / 4 W y/25 



where B(l) is standard Gaussian. As remarked by Arratia [3, p. 71], the conclusion for general 
5 > follows from the case 25 = 1 by Brownian scaling. See also De Masi, Ferrari |9], Rost, Vares 



|39| and Arratia |3] where variants (or generalizations) of (1.2) are proved for a tagged particle in 
the exclusion process on Z associated with a simple symmetric random walk. Harris conjectured 
that the process Br ) is not Markov, and left open the problem of describing the long-run behavior 
of the paths of Br y As far as we know, these problems are still open. 

The case of distinct drift parameters differs from these classical models in several remarkable 
ways. For example, if I = {1, . . . , K} for some K £N and 

1 j 1 K 

J k=l J k=j+l 

for all j = 1 , . . . , K — 1 , then there exists an invariant distribution for the process of gaps between 
consecutive particles (X^ 2 )(t) — Xn\(t), . . . ,X^\(t) — X(K-u(t)), t > (see Theorem 8 in |37j ) 
which is not the case if 5± = ■ ■ ■ = 5k- When I = N, one can consider the so-called Atlas model in 
which 5± = 5 > and 5i = for every i ^ 1 in I. In this case it is shown in |37| that, if the initial 
positions of the particles are chosen according to a standard Poisson process of rate 25, then the 
joint distribution of the gaps between consecutive particles Xf 2 )(t) — Xn-jit), Xr s \(t) — X^(t), . . . 
is the same for all t > 0. However, the same statement is not true when 5{ = 5 for all £ € J (see 
Theorem 4.2 in |40|), Moreover, if we consider the motion of the left-most particle X^(t), t > in 
the Atlas model, no estimates on its growth and fluctuations are known. Our article is a step in the 
latter direction. Using techniques from the theory of concentration of measures we give estimates 
on fluctuations of the paths of the distances between ordered particles and associated statistics. 

One such statistic is given by the market weights. The latter can be defined for I = {1, . . . , K} 
with an arbitrary K G N and any choice of 5\, ... ,5k by setting 

e X( K _ i+1) {t) 

Hi(t)= K {t y i = l,...,K (1.3) 

for all t > 0. It is clear that for any fixed t the sequence fii(t), . . . , f^K(t) is a non-increasing 
sequence of positive numbers that add up to one. These numbers, in econometric models, go by 
the name of market weights and have an interesting history. 
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Figure 1: Capital distribution curves: 1929-1999 



Fernholz in his 2002 book [19] introduces solutions of (1.1) to model the time dynamics of the 
logarithmic market capitalizations of different companies in an equity market. In other words, he 
considers a stock market with K companies whose total worths in stocks are given by exponentials 



of the one-dimensional components of the solution of equation (1.1 ). A major objective of his work 
is to explain the following curious empirical fact. Consider the shape one obtains by plotting log fa 
versus logi. This log-log plot is referred to as the capital distribution curve. See Figure 1 below 
(reproduced from |19| ) which shows the capital distribution curves between 1929 and 1999 for all 
the major US stock markets (NYSE, AMEX and NASDAQ) combined. Empirically, the left part 
of the plot exhibits nearly linear decay. This corresponds to the market weights, in decreasing 
order, displaying power law (or Zipf 's law) decay. More strikingly, the slope of the decay is nearly 
constant over eight decades, something truly remarkable in the volatile world of financial markets. 

In [7] the authors explain the linearity by proving (under suitable assumptions) that the possible 
limiting stationary laws of the market weights, as K grows to infinity, are given by a subset of the 
Poisson-Dirichlet family of distributions. The masses of this family have a polynomial decay with 
size, which corresponds to the linear decay in the log-log plot. However, this does not quite address 
the temporal stability of the picture above. Our analysis below captures some of its subtleties. 

To analyze the stability of the shape of the capital distribution curve, consider the process of 
market weights (fa(t), i = 1, . . . ,K), t > and fix a J <C K. At any point of time t > we 
introduce the linear regression between the pairs of data {(log i, log fa (t)), i = 1, . . . , J}, passing 
through the first point (0, log^i(f)). In other words, we deal with regressions of the form 

log fa(t) -log fa(t) = a(t) log i + ei{t), i = 2,...,J (1.4) 

for fixed values of t > 0. Clearly, the resulting ordinary least squares estimator for the slope 
parameter corresponds to the slope of the linear part of the curve at time t in Figure 1. With a 
minor abuse of notation we will denote this estimator by a(t). 

Our objective is to estimate the fluctuations of the above parameter when the spacing process is 
running close to its stationary law. We choose the Atlas model to have a specific sequence of drifts 
8%, ... ,8k, although other values of the drift parameters can be easily substituted. We consider 
initial configurations the spacings of which are close to their unique invariant distribution found in 



|37j (see Lemma 3.1 below). 
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Theorem 1.1 Let I = {1, . . . , K} and consider the Atlas model in which 5\ = 5 and 5i = for all 
i = 2, . . . , K. Define the initial spacings between the particles by 

X {K) (Q)-X {K _ i+l) (0) = jlogi, i = 2,...,K. (1.5) 
Assume that K is sufficiently large. Fix a J < K/15 and consider the process a{t), t > of 



ordinary least squares estimators resulting from regressions in (1.4) for different values oft. Then 
the process a(t), < t < 5~ 2 K satisfies the following concentration of measure property. Let 
a = sup < s < s -2 K [a(s)}. 

Then there is a constant m n such that 



F(a<m a )< l/2 + nVKe- K/500 , 
and for all positive r sufficiently large, one has 

F[a> m a + rVKj < 2exp . 

Hereby, fi is an absolute positive constant and C a is a positive constant depending on J and given 
by 

'EffiogC-zi/ioV 



C a (J) = J 



3 



Et 2 (iogi) 2 



Remark. The above result, although novel, does not capture fully the extent of concentration that 
is seen in the real world data (see Figure 1). We suspect that the reason for this is the empirically 
observed unequal diffusion coefficients of the ordered particles. See, for example, the discussion 
in [3] which mentions that the diffusion coefficient for the ranked particles decays linearly with 
increasing rank. Thus, particles at the top are more stable than the average which should lead to 
a higher concentration. 

Remark. Another obvious way to improve the bound is to use true stationarity. One can divide 
up a long time interval in small subintervals of appropriate size, use the above theorem on the 
small intervals, and take a union bound. However, in this case it is not clear if the Gaussian 
concentration is preserved when the invariant distribution has Exponential tails. However, we 
expect that variations of our method and the argument in this remark can lead to Exponential 
concentration over much larger intervals. 



Next, we set / = N in (1.1 ) and consider a sequence of drift coefficients 61,62, ■■ ■ which satisfies 

6m = $m+i = • • • (1-6) 

for some M G N. Under the assumption that the sequence of initial positions of the particles 
X\ (0),X2(0), ... is deterministic, non-decreasing and such that 

liminf^Mx), (1.7) 

i— >oo i 



the system of stochastic differential equations (1.1) has a unique weak solution (see Proposition 3.1 



in |42j). We let £(1^+1) (i), t > be the local time process at zero of the process Xu+i\(t) — X^(t), 
t > for each i G N and are interested in the concentration properties of the vector-valued processes 
(I (2) (t)-I (1) (i),..,I w (t)-I M) (t)),te[0 ) T] and (L {1>2) (t),...,L [n _ 1>n) (t)),te [0,7] for 
arbitrary values of n G N and T > 0, which we consider to be fixed from now on. For our main 
result on the latter we introduce the following assumption on the initial particle configuration 

(x l (o),x 2 (o),...y. 
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Assumption 1.2 The sequence of initial positions of the particles is non- decreasing with probability 
one, and there exists a deterministic constant c > such that 



X k (0)-X N (0) >c{k-N) 
holds for all k > N almost surely, where N = max(n, M). 



a. 



We note that Assumption 1.2 ensures the existence of a unique weak solution for the system 
fllT] ) by Proposition 3.1 in 02]. 

To state our second main theorem we define the norm 

n-1 



T,2 



1 



n 



sup fi{t) 
1 7~r o<t<r 

^=l — — 



l 

2\ 2 



(1.9) 



on the space C([0, T], M n of continuous W 1 1 -valued functions on [0, T] where fa, i = 1, . . . , n— 1 
are the component functions of /. Our result then reads as follows. 



Theorem 1.3 Let the Assumption \1.£\ be satisfied with a constant c > 0. Moreover, let A and B 
be measurable subsets of (C([0, T],R n ~ ), ||.||t,2) suc/i that 



P(((L (M +i)(t),» = l,...,n-l), f£[0,T])Gi) > i 

i+1) (t)-X (i) (t),i = l,...,n-l), tG[0,T])£5) > i 



and /or any r > set 



Ar = {fteC([0,T], 
B r = {he C([0,T], 



t>n—l\ 



r>n—l\ 



jnf || h — h\\T,2 < r } ; 

jnf ||/i — h\\T,2 < r }- 



(1.10) 
(I.ll) 

(1.12) 
(1.13) 



T/ien t/iere exist positive constants C\, Ci, C3, C4 depending on c, A = max,=i ... jvf-i \5j — $m\, 
M, n, T and the value of the left-hand side of (1.10) and (I.ll), respectively, such that for all 
r\ > C\ and r^ > C2 it holds 

I (n — 1)7 c 7 



F[{(L ii;i+1) (t),i = l,...,n-l),te [0,T])<£A ri ) <C 3 exp 
F(((X (i+1) (t)-X {t) (t),i = l,...,n-l),tG[0,T])^B r2 ) <C 4 exp 



14T 

I (n — l)7c 7 
" r2 257^ 



Remark. The local time of collisions between two consecutive ordered particles, as considered in 



Theorem 1.3, is interesting both mathematically and in applications. Mathematically say in the 
Atlas model, the local time (compared to the the case Si = 0, i € I) measures the push felt by the 
various particles due to the drift at the bottom. Its significance in economic models is discussed 
in Section 13 of [20J, which also mentions the somewhat surprising fact that these local times can 
indeed be measured from data. 

The rest of the article is organized as follows. In the next section we recall some facts about the 
concentration of measure phenomenon and the Skorokhod problem in the orthant, and provide an 
upper bound on the Lipschitz constant for the Skorokhod map of interest in Lemma [2. 4| The latter 
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is the key to the proofs of the two main results. Its proof relies on the construction of the Skorokhod 
map by Harrison and Reiman |24j and applies to other Skorokhod problems in the orthant as well. 
In Section [3] we use a relation between Transportation Cost Inequalities and the concentration of 
measure phenomenon to complete the proof of Theorem |1.1| and provide the remainder of the proof 
of Theorem 11.31 in Section HI 

2 Preliminaries 

2.1 Some facts about measure concentration for stochastic processes 

Suppose (X , d) is a complete separable metric space equipped with the Borel cr-algebra. For a Borel 
subset A of X and a positive real number r, define the r-neighborhood of A by 

A r :={x£X: d(x,A)<r}. 

We say that a probability measure [x on (X,d) has the measure concentration property if for any 
Borel subset A with [J>(A) > 1/2 the value of n(A r ) is very close to one for large values of r. This 
closeness is usually expressed as a sub-gaussian tail in terms of r. 

The concentration of measure phenomenon has become one of the most important concepts 
in modern probability theory. For an excellent introduction to this area we refer the reader to 
the beautiful article by Talagrand |47j . A considerable effort has been spent by probabilists on 
identifying distributions that have the measure concentration property. A (somewhat dated) survey 
can be found in the monograph by Ledoux |29j . 

One technique for proving the measure concentration property, originally proposed by Marton 
|30| . involves the so-called Transportation Cost Inequalities (TCI) that we describe below (see also 
Talagrand [3S], Dembo |10j . Dembo and Zeitouni |llj). Consider, as before, a complete separable 
metric space (X, d) endowed with its Borel cr-algebra. For a real number p > 1 and all probability 
measures P and Q on the latter space define the p-th Wasserstein distance 

W P (P,Q)=M [Ed(X,X') p ] 1/p , 

where the infimum is taken over all couplings of a pair of random elements (X, X') such that the 
marginal law of X is P and that of X' is Q. 

Next, fix a particular probability measure P. Suppose there is a constant C > such that for 
all probability measures Q<Pwe have 

m(P, Q) < V2CH (Q | P) (11.14) 

where H refers to the relative entropy H (Q | P) = log (dQ/dP). In this case we say that P 
satisfies the Quadratic Transportation Cost Inequality (QTCI) with the constant C. 

A function / : X — > K will be called Lipschitz if there is a positive constant a for which 

1/0*0 - f{y)\ < ad(x,y), x,yeX. 

The smallest such constant a is then referred to as the Lipschitz constant of /. We shall call a 
function 1-Lipschitz if a can be taken to be one. Let C denote the set of all 1-Lipschitz functions 
on (X, d). The (very short) proof of the following theorem can be found in Ledoux |29[ p. 118] and 
the original article by Marton |30| . 

Theorem 2.1 Suppose that P satisfies the QTCI with constant C. Then one has the following 
concentration estimates for all r > 2\/2Clog2: 
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(i) For any measurable set A such that n(A) > 1/2, it holds 



H(A r ) > 1 - exp (-r 2 /8C) . (11.15) 

(ii) For any f 6 C, one has 

P{x: \f{x) -m f \>r) <2e~ r2/8C , (11.16) 

where mj is the median of f with respect to P. 

The QTCI is unique in its advantages and is related to the log-Sobolev inequality, hypercon- 
tractivity, Poincare inequality, inf-convolution and Hamilton- Jacobi equations. For details, please 
consult Otto and Vilani |35| . Bobkov and Gotze [5], and Bobkov, Gentil and Ledoux [6j. 



In addition to Theorem 2.1 the following lemma will be useful in the later text. Its (short) 



proof can be found in Lemma 2.1 of |12j . 

Lemma 2.2 Suppose /x is a measure on a metric space (E,oIe) that satisfies the QTCI with a 
constant C. Let (F,d F ) be another metric space. If the map ^> : (E,dE) — > (F,d F ) is Lipschitz, 
i.e. 

d F {^{x),^{y)) < ad E (x,y), x,y € E 
for some constant a > 0, then p, = fj, o satisfies the QTCI on (F, d F ) with the constant Ca 2 . 

In this article we are interested in the choice of X = C([0,T],M i ^- 1 ), the space of continuous 
maps from the interval [0, T] to where K is the number of interacting particles as before. The 

latter space is typically referred to as the path space. It will be endowed with different variants of the 
uniform metric which are described below. We shall use Quadratic Transportation Cost Inequalities 
satisfied by the laws of stochastic processes with continuous paths, especially multidimensional 
reflected Brownian motions. 

Let us provide a brief review of the literature on measure concentration in path space. Houdre 
and Privault |25j . and Nourdin and Viens |34| use tools from Malliavin Calculus to derive concen- 
tration inequalities for functionals on the Wiener space, among other things. Transportation Cost 
Inequalities on the Wiener space first appeared in Feyel and Ustiinel [21]. Djellout, Guillin and Wu 
|12| provide a characterization of the L^TCI for diffusions. They also prove the QTCI for diffusions 
with respect to the Cameron-Martin L 2 -metric. Several articles in analysis and geometry are also 
devoted to this topic: Fang and Shao [HI H7], Fang, Wang, and Wu |18j . Gourcy and Wu |22| . 
Wang |50| I51j . and Wu and Zhang |52j . The QTCI for diffusion-like semimartingales and reflected 
processes with respect to the uniform distance appeared in Pal |36| and was generalized by Ustiinel 
|49j . We refer the reader to |36| for more details on the existing literature. Similar ideas have been 
also put forward for discrete Markov chains. See, for example, the articles by Marton [31[ 1 



2.2 RBM and the Skorokhod map 

The reflected Brownian motions (RBMs) we are interested in have a constant drift vector, a constant 
diffusion matrix and are reflected whenever they hit the boundary of the positive orthant. On each 
face of the boundary of the latter the direction of reflection is constant. The theory of such processes 
is well-developed. In particular, their existence and pathwise uniqueness follows from the existence 
of a deterministic transformation mapping Brownian paths to the corresponding reflected paths. 
This is the so-called Skorokhod map, whose one-dimensional version is due to Skorokhod. We 
lift the following description from the article by Harrison and Reiman |24j . For more details and 
generalizations see [13], [2], and |15| . 
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Define C([0, oo), l^ 1 ) as the space of continuous functions x : [0, oo) — > ~$L K 1 , endowed with 
the topology of the locally uniform convergence. For each such function x we denote its component 
functions by x% for i = 1, ... ,K — 1. Let Cs denote the subset of functions in C([0, oo), 
for which Xj(0) > 0, i = 1, ... ,K — 1, and let Q = {qij) be a (K — 1) x (K — 1) matrix with 
non-negative entries, zeros on the diagonal and spectral radius strictly less than one. 

Theorem 2.3 Q24J) For each x G Cs there exists a unique pair of functions y,z 6 C([0, oo), R^ -1 ) 
satisfying 



K-l 



8=1 



zAt) > 0, t > 



(11.17) 
(11.18) 



for all j = 1, . . . , K — 1 an<i suc/i t/iat for every i = 1, . . . , K — 1, the function yi increases only at 

those times t for which Zi(t) = 0. 

Moreover, suppose that the matrix Q satisfies 



K-l 



WQWcs '■= . , Yl qi 3 < L 

j=l,...,K-l ^ 



1=1 



(11.19) 



Then y is the unique function in C([0, oo), M^ -1 ) that for all t > satisfies the vector equation 
y(t) = sup < s <£ [y{s)Q — x{s)] + , and is given by the limit (in the locally uniform topology) of the 
following iterative scheme: 



y [0 Kt) - o, 

y^ i+1 \t) = sup ly m {s)Q-x(s) 



t > 0. 



0<s<t 

Hereby, the supremum and the positive part are taken componentwise. 



(11.20) 
(11.21) 



It will become apparent later that the matrix Q of our choice will be given by 



-i) 



/ o 

1 

2 







V o 



i 

2 





o \ 





1 

2 

o / 



Since Q^" 1 ) is a finite, irreducible, substochastic matrix, it is immediate (by adding an ab- 
sorbing point) that the spectral norm of Q^ K ~^ is strictly less than one. As was shown in |24] , 
this implies that the corresponding Skorokhod map, which transforms the input path x into the 
reflected path z (or, the path of the "local time" y), is Lipschitz with respect to a variant of the 
\>T2 defined in (1.9). We would like to obtain an explicit upper bound on its Lipschitz 



norm 



constant with respect to the norm ||.||t,2- 

Note that Q^ K ~ l "> does not satisfy the Assumption (11.19) and, hence, the iterative scheme of 



Theorem |2.3| cannot be applied directly to construct the corresponding Skorokhod map. The way 
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to get around this is to define a (K — 1) x [K — 1) diagonal matrix 



£>(*-!) 



(K-l) 2_ 
K 



\ 



w (K-D(K=k) J 



with a strictly concave function w^ K ^ : [0,1] — >• [0,oo) such that w^ K ^(u) = if and only if 
u G {0, 1}. Viewing the set of equations in ( |II,17[ ) for different values of j as an equation between 
row vectors, and multiplying both sides of it by the matrix [D^ -1 )] from the right, we obtain 
the equations 



K-l 



z ' 3 {t) = x'^t) + y'^t) - njyKt), t > o 



(11.22) 



8=1 



for j = 1, . . . ,K - 1 where x' = x [D {R - 1 )] 1 , y' = y [D^-^] \ 

R^-D = (r ..) = D (K-1) Q (K-1) l D (K-l 



z = z 

) - 1 



[ D (K-1) 



-1 



and 



(11.23) 



Note that the coordinates of x' , y' and z' can be computed by a simple rescaling of the corre- 
sponding coordinates of x, y and z, respectively. Moreover, it holds 



w 



(K-l) (1=1 



max 
1=1,..., K-l 



+ w (K-i)(m 



2w^-i)U 



< 1 



(11.24) 



due to the strict concavity of the function w^ K l '. Thus, the Skorokhod map corresponding to the 



matrix 

R(K-1) 

can be obtained using the iterative scheme of Theorem 2.3 



From now on we fix a terminal time T > and equip the space C^T],!^- 1 ) of continuous 
^ _1 -valued functions on [0,T] with the norm 



/ 1 K-l 

\x\\t,2 ■= 7 ^2 SUp Xi(t) 

\ K ~ 1 ~[ 0<t<T 



(11.25) 



where X{, i = 1, ... ,K — 1 are the component functions of x. We write & L and <I>^ 

for the maps that take a path x G C([0, T], l^ 1 ) to the local time path y G C([0, T], M^" 1 ) 
and the reflected path z G C([0, T], M^ -1 ), respectively, corresponding to the Skorokhod problem 
with reflection matrix Q( K ~^ defined above. The following lemma is one of the crucial steps in 
the proofs of Theorems |1.1| and 1.3 It provides an upper bound on the Lipschitz constant of the 
Skorokhod map and is a significant improvement on an earlier attempt by Pal |36j . 

Lemma 2.4 For all natural numbers K > 2 the map 

$2 is Lipschitz on C([0, T], M^ 1 ) with 



respect to the norm ||.||r,2 defined in (11.25). Moreover, its Lipschitz constant Lip L 

Lipf- l) < 2- (if - 1)5. 



(K-l) 



satisfies 
(11.26) 



9 



Proof. 1) We fix a natural number K as in the statement of the lemma and will prove the inequality 
(11.26) for that value of K. From the considerations preceeding the lemma we know that for each 
g G C([0, T], R^ _1 ) the map K ^ corresponding to the Skorokhod problem with refle ction 

with 



2.3 



matrix R^-V, evaluated at g, is given by the limit of the iterative scheme of Theorem 
input g, restricted to the interval [0, T]. Now, let g be another function in C([0, T], R^ -1 ), for each 



k £ N define f k and f k as the results of the k-th step of the iterative scheme of Theorem 2.3| witrj 
inputs g and g, respectively, and set / and / for K 1 (g) and ^¥~ K ^ (Jf), respectively. Finally, 
define the norm ||.||T,max on C([0, X 1 ], IR.^" ^) by 

\\x\\T,max= max sup \Xi(t)\, (11.27) 

i=l,...,K-l o<t<T 

where Xj, i = 1, . . . , K — 1 denote the component functions of a function x 6 C([0, T],R —1 ) as 
before. By applying the triangle inequality, the fact that the operation of taking the positive part 
is 1-Lipschitz and the definition of the norm \\.\\ C s in (11.19), one obtains the following chain of 
inequalities: 



with 



T,max 



max sup 

i=l,...,A'-l o<t<T 



sup ((f k (s)R ( - K -% - g t (s)) + - sup ((f k (s)R^ K -% -&(*))- 

0<s<t 0<s<t 



< max sup sup 

i=l,...,K-l Q<t<T0<s<t 



((f k (s)R( K -%-g t (s)) + - ((JH^-^-U^ 



<. max ( sup \(f k (t)R( K -%-(f k (t)R( K -%\+ sup \g t (t) -g t (t)\) 

t=l,...,K~l \o<t<T 0<t<T 7 

^H^-^llc-. max sup |/f(t)-^(t)| + || 5 -5|| Timax . 
i=l,...,K-l o<t<T 

Taking the limit k — > oo and rearranging terms we conclude 

11/ - f\\T,max < 2 _ ' ' I' 5 ~~ 9\\T,max- (11.28) 

Recalling that the map $^ can be obtained from the map <3?^ by a rescaling of the 
coordinates of according to the matrix D^ K ~ 1 ^ (see equation (11.22)), we deduce 

W^T^id) - $f K " 13 \g)\\T,max < !_ II^-DH^ • \\9-9\\T,max, (H.29) 



(11.30) 



where 

maxi = i i ... )i <r_iu;( K - 1 )( : ^) 

Q = t — ■ 

min /=li ... jA '_ 1 w( x - 1 )(^) ' 

Thus, the obvious equivalence of norms inequalities between the norms II. II y max ^rid II. II 7*2 on 
C({0,T\,'R K - 1 ) show 

Q 11 ~,| ^ (^-l)^g I, ~,| 

- HpP|i " 115 " 5l|T ' mra - i-||^-D|u ' 115 " 5||t ' 2 " 

Since the functions g and were chosen arbitrarily in C([0,T],M^ _1 ), we conclude 

Lwf^ < ^'^Z • (H.31) 
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2) For K = 2 or K = 3 one can take D^ K ^ to be the identity matrix of appropriate dimen sion, 



so that q = 1 and \\R^ K ~ l 



\\Q^ K ^'\\cs which immediately yields the inequality (11.26). To 



finish the proof for the case K > 4, we need to bound the right-hand side in (11.31 ) from above for 
a suitable function w( K ~ l > in the definition of the matrix D^ -1 ). From (11.24) it is clear that 



(x) + w 



(K-l) 



< 



max 
0<s<l-- 



x + 



K 



2w (K-l)( x+ j_ 



(11.32) 



To define w^ K ~ l \ we choose a continuous function v^ K ~ 1 ' ) on [0,1] such that v^ K ~^{^) = 0, 
v( K ~^(x) = — v^^'^l — x), x £ [0, 1] and t/^ -1 ) is strictly decreasing and convex on [0, h]. Then 
we set w^^^fyx) = Jq v^ K ~ 1 \y) dy for all x S [0, 1]. In this case the convexity properties of 
and w( K ~ ^ imply that 



d w^ K - l \x)+w { - K -^[x + j ? ) 
dx 



2w (K-l)( x+ j_ 



> 



(11.33) 



for x £ [0, K — j(\ and that the opposite inequality holds for x G [g, 1 — jA- Moreover, if the latter 
derivative was zero at a point i£ [k — j^, I — -^), then the strict concavity of the function 

would imply i + ^ ) < 2u^ _1 Ux + — f^ _1 ^(x). However, this is a contradiction to 



the fact that the graph of the function v^ K ^ on [i, I + is never below the line connecting the 



points (giO) and 



+ 



K ' 



+ 



A' 



All in all, we conclude that the maximum in 



w 



K 



max 
l=l,...,K-l 



is achieved when 4= takes the value closest to \. If if is even, the value \ is attained and, thus, 



the value of the latter maximum is given by 

+ w (K-l)(l + l 



It 



(K-l) 1 _ J, 
2 K 



K 



2wi K - 1 )(\ 



1 



,(K-1) 



2 K 
1 



(y) dy 



gv^Kv) dy 

If K is odd, the maximum is achieved at 



jh k v (K-l) {y) dy + k 



5 + ..(K-l) 



(y) dy 



2jiv^){y)dy 



^yJ 1 and for its value one has 



< i 



Ji^v^Xy) dy 

2K 

2j}v(K-i)(y) dy' 



2K 



IX — 1 I\ -f- 1 

J3 v^Hy) dy - f& t^Jfc,) dy 

2 A' 'Hi 

2 f ^r v (K-V(y) dy 
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Hereby, the last inequality is a consequence of the monotonicity and the symmetry properties of the 
function v^- K ~ l \ which show that we have decreased the numerator and increased the denominator 
of the fraction. The latter two calculations show that 

l \ 2 



-i uJk-du ^ T~i \ — / — 1 — " V (IL34) 

1 " I|jR( )||cs v(K-V(y) dy) ■ ( J* vi*-i)(y) dy) 

For the unique function v^ K ~^ of the type described above for which v^ _1 ^(0) = K and v^ K ~^(^) = 
2, and which is affine on [0, j^] and [j^, the latter expression computes to 

<2-(K-l) 2 . (11.35) 



2 K+2 
(K~2)K ' 2K 



We remark at this point that the order of magnitude of the bound in terms of K is optimal for 
this choice of the function v^ K ~^\ since the non-negativity and the convexity of v^ K ~ 1 ^ on [0, \] 
together with i/^ -1 )^) = show 

K v^ K ~ l \y) dy < JJ v^(y) dy, j" ^ v^^y) dy < {^f ■ v^Xy) dy. 



2 K 



Finally combining the inequalities (11.31), (11.34) and (11.35), we end up with the statement of the 



lemma. □ 

3 Concentration of the shape of the market weights 



Consider the particle system in (JTTTj) with / = { 1 , . . . , K } for some if £ N and a choice of constants 
Si,..., 6k that satisfy the following condition. Setting 6 = -4 Y2i=i S%> one nas 

j 

«i : = " S K-i+i) > 0, j = 1, . . . , K - 1. (111.36) 

i=l 

In our analysis of the market weights we will use the following result from |37j . 



Lemma 3.1 Under the condition (IIII.36 ) the process of spacings 



(£i(i): i = l,...,K-l):=(X iK _ l+1) (t)-X {K _ i) (t): i = 1, . . . , K - 1), t>0 (111.37) 

has a unique stationary distribution which is that of independent exponential random variables with 
rates 2a i, for i = 1, . . . ,K — 1. Moreover, the system of spacings is reversible at equilibrium. 

A situation in which the above lemma applies is given by the Atlas model with 6i = 5 > and 
6i = for alH = 2, . . . , K. In that case one easily computes 

* = i = l,->*-l- (HI.38) 



i=l 



Thus, under the stationary distribution we have 

K 
26~j- 

j—l 

nX(K-i+i)(t) - x [K _ j+1) (t)} = ^ ^ y « - log(iA), i < * < i < # 
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for all t > 0. 



Next, consider the linear regression in (1.4 ). The corresponding ordinary least squares estimator 
for the slope parameter is given by the formula 

_ Ef= 2 (l°g*) (log/Ji(t)-log/Xi(t)) _ Ei=2^Si)(X (K) (t) - X (K _ i+1) (t)) 

Ei=2 lo g 1 Ei=2 lo g 1 

E I J =2 ( 1 °gO ggkjg) = Eg^gtj+ijogj fIII 39) 

E/= 2 iog 2 ^ Ei 2 iog 2 ^ 
_ E/=i 1 Qog^A0^) 

E£V(i + l) 
We use this formula in the proof below. 

Proof of Theorem \l.l\ The proof is broken down into several steps. 

Step 1. A Quadratic Transportation Cost Inequality for the process of spacings is derived in the Atlas 
model with K particles. 

Step 2. Next, we assume K to be very large compared to J. We prove a localization lemma that 
shows that the process a(t), < t < 5~ 2 K is determined only by the particles corresponding 
to the 5J topmost indices with very high probability. 

Step 3. Finally, we show that the law of the 5 J particles in Step 2 is approximately that of another 
rank-based process, so that we can use the estimates obtained in step 1, with K replaced by 
5 J, to bound concentration of measure probabilities under the event constructed in step 2. 



Step 1. Consider the process • ■ ■ , Xrj(\(t)), t £ [0, T] of ordered particles in the system 

(1.1 ) with / = {1, . . . , K}. From Lemma 4 in [37] we know that there exist i.i.d. standard Brownian 
motions Pi, ... , /3 k such that for alii = 1, . . . ,K it holds 

dX {l) (t) = Si dt + d/3i(t) + -dL {i _ u) {t) - -dL (i)i+1) (t), t e [0,T] (111.40) 

with Li t i + i(t), t G [0, T] being the local time process at zero of the process Xr i+ i\(t) — X^(t), 
t G [0, T] for i = 1, . . . , K — 1 and the convention Lo,i(i) = Lj^k+i (i) = for all t E [0, T\. 
Hence, for all i = 1, . . . , K — 1 one has the dynamics 

d{X {l+1) {t) - X {i) (t)) = (S i+1 - Si)dt + df3 i+1 (t) - dPiit) 

+ dL(i,i+i)(t) - -dL(i +1}i+ 2)(t) - 

on [0, T}. In particular, we can conclude from this representation as in section 2 of [37J that the 
process 

(X {2) (t)-X il) (t),...,X {K) (t)-X {K _ 1) (t)), te[0,T] 

is a reflected Brownian motion in the (K — l)-dimensional positive orthant with reflection matrix 
q(^-!) i n the sense of section 1 in [24] . By reversing the labeling we see that the process := 
. . . , ^K-i(t)), t G [0, T] is also an RBM in the positive orthant with reflection matrix Q( K—1 \ 



13 



Thus, the process £ can be obtained as the image of the process 7*(t) := + S^it — f3i(t) — 

Sit, i = 1, . . . , K - 1), t G [0, T] under the map K_1> . 

By Theorem 6 in [36j the process of independent Brownian motions (Pi(t)+5\t, . . . , /3jc(i)+SKt), 
t G [0, T] satisfies a QTCI with respect to the norm ||.||t2 with the constant 4K~ 1 T. Moreover, 
the map that takes the vector of these Brownian motions to the process J*(t), t G [0, T] is L ipsch itz 



, 2 with Lipschitz constant 2yJ Kj (K — 1). Hence, by Lemma 

with the constant 



with respect to the norm 

process J*(t), t G [0, T] satisfies a QTCI with respect to the norm ||-|| r2 

K 



2.2 



the 



C* K (0) ■= WK^T 



K-l 



16(K- l) -1 ^ 



(111.41) 



Now, we will use Lemma 



2.4 



Consider the equation (11.17) with the matrix Q( k ^ for a fixed 

>(if-i), \ 



value of t and two different unconstrained processes x and x. Writing y, y, z, z for & L 
<I>^ 1 \x), (x), > (^)) respectively, we easily deduce 

K-l 

\zj(t) - Zj(t)\ < \ Xj (f)-Xj(t)\ + ( . max |y< - ^| ) • V I (J^" 1 ) - Q^" 1 ^^- , tG[0,T] 



i=l 



<2 



where x ) is a (X — 1) x (ET — 1) identity matrix. Hence, 

lT,max — IF — ^llT,max + ^\\V ~ mT,n 



Z — Z 



On the other hand, combing the bounds (11.29), (11.34) and (11.35) we get 



(111.42) 



|y-y|lT,max<2(^-l) 2 ||x 



X\ 



T.max 



(111.43) 



for all K > 4. By the same argument as in the beginning of step 2 in the proof of Lemma 2.4 the 
same inequality is true for K = 2 and K = 3. Putting the inequalities (III. 42) and (III. 43) together 
we obtain 



z — z 



T,max 



< (1 + 4(K-1) 2 ) 



x — X 



T.max 



< 5(K-1) 



x — x 



T,max 



(111.44) 



for all K > 2. Recall that by applying the map to the paths of the process 7* and by 

reversing the order of the coordinates thereafter one gets the process £. Combining this observation 
with Lemma 2.2, (III. 41) and (III. 44), we see that the process £(i), t G [0,T] satisfies a QTCI with 
respect to the norm ||-||^ max with the constant 

C* K := 400(if - 1) 3 T. (111.45) 

Now, we restrict ourselves to the first (J — 1) coordinates of the process £, since only those 
appear in (III. 39). By Lemma 2.2 the vector-valued process (£i(i), • ■ ■ , £j_i(t)), t G [0,T] also 
satisfies a QTCI with respect to the norm ||"||y max with the constant C* K . 



Step 2. Consider the formula for a(t) in (111.39). The value of a(t) depends only on the top J 
spacings, whereby we have assumed that J is very much smaller than K. In this case, for a large 
enough m and with high probability, the top J processes during the time interval [0, T] are identical 
to the top J processes among the processes which start off at the top J + m positions at time zero. 
The following lemma makes this idea precise. 
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Lemma 3.2 Consider the particle system of Theorem 1.1 and for all m 6 N define a m as the first 
time t at which, for some i > J + m and some 1 < j < J, it holds Xx-i+i(t) = X^-j+i){t)- 
Then, for all K > 30, we have 

P (cj4j + i < S- 2 K) < 6J 2 K- 3/2 e~ K/50 . (111.46) 

Proof of the Lemma. For any T > 0, the event {o~ m < T} implies that for some i > J + m and 
some 1 < j i < J, the processes and X^-j+i cross paths during the time interval [0,T]. 

Using the union bound and bounding the drift of the lower particle by the constant 5, we have the 
following estimate: 

P (c m <T)< V V P f sup (W K -i+i(t) - W K -j+i(t)) >-5T + X K _ J+1 (0) - XK-i+ttp)) 
J V P ( sup (W K . i+1 (t) - Wk-j+iH)) > -5T + X K „ J+1 (0) - X K -i+i(P)\ • 

i=J+m \ / 

(111.47) 

Next, we use the fact that the supremum of a standard Brownian motion up to time T has the 
same law as the absolute value of a normal random variable with mean zero and variance T. Thus, 
with a standard normal random variable Z one has 

P( sup (WK-i+i(t)-W K -j+i(t))>-6T + X K „j +1 (0)-X K - i+ i(0) ) 

\0<t<T J 

/2f\Z\>j log(VJ) - STj = 2$ ^ \og(i/J) - S^j j . 

Here, $ is the one minus the cumulative distribution function of a standard normal random 
variable. 



< 



Plugging this into (III. 47 1, we get 



>(a m <T)<2J f; ^((iW)-^) ) 

i=J+m \ \ V / +/ 



oo 



< 2J I *l [j^^ g (x/J)-5^) )dx. 



J+m-l 



(111.48) 



+/ 

We note that for m > 4 J + 1 and T = 5~ 2 K we may omit taking the positive part in the latter 
formula. 

Now, recall the well-known inequalities (see [TJ p. 298]) 

2lf{ jj— <Hy) < -<p(y), y>0, (111.49) 

where (f is the density function of a standard normal random variable. From the latter inequality 
it follows that the function log$(y) is concave on (0, oo). To see this, note that 



15 



In other words, for any positive real numbers a and b, we have 

-,- ,'a + b" 



> */$(a)$(&). 



Taking T = 5 2 A and m > 4 J + 1, and choosing 



if /T /T 

" = _ ; — \og(x/J) - 8\ -, b = Sy 



6V2T V 2 

for x > J + m — 1 in the latter inequality, we get 



* & og(I/J) - s vlp 



A" 



25V2T 



log(x/J) 



for all x > J + m — 1. Substituting in (|III.48|), we deduce 



(o- m < T) < 2J 

4J 2 (K/2T 



n -1 



A' 



J+m-l 
1 



25V2T 



log(ac/J) 



r/.r 



exp [2SV2fy/K) $(y) 



where a* = log[(J + m — 1)/J] K/(2Sy2T). Hence, using < y 1 f(y), y > 0, we obtain 

-l 



1 (Om < T) < 



< 



2J 2 5V2T 
irKa* 2 

2J 2 <5 V / 2T 
irKa* 2 



$ Si 



$ S\ 



2 



a* 



exp (25V2Ty/K -y 2 ) dy 



e 2Tsy K ^ a* 



5V2T 
A 



At this point we take m = 4J+1 and substitute 8 2 K for T. In this case, one has a* 
and the latter bound simplifies to 



{(TAJ+l < S- 2 K) < ( 1 



A 3 / 2 



-1 



Hereby, we have set 



16^2J 2 2K -i 



7r(log 5) 2 



Using the bounds in (III. 49) we obtain 



2A-V2 



1/V2+ v/V2 + 4A- 1 



teVK- -2= 



Coy/ A 



16 



for all K > 30, where 

~ log5 2 

Hence, combining all our estimates, we end up with 



p k, +1 < i-A-) < ^ . j — exp y~4 ~ ~~2 j ' 

To finish the proof of the lemma, it suffices to observe that the right-hand side of the latter inequality 
is bounded above by 6J 2 K~ 3 / 2 e- K / 50 for all K > 30. 

Step 3. Consider the particle system of Theorem |1.1| We proceed with another localization lemma. 

Lemma 3.3 Define a m to be the first time t at which, for some i > m, we have Xi(t) = Xr\\{t). 
Then 

F(a 2K/3+1 <6- 2 K)<llVKe- K / 500 . 

Proof of the Lemma. The proof is very similar to that of Lemma 3.2 and we only outline the 
argument. With m = 2K/3 + 1, and letting Z be a standard normal random variable and $ be 
one minus its cumulative distribution function as before, we have 

F(a m < 5- 2 K) < VP [ sup (Wi(t) - Wi(t)) > -S~ X K + X t (0) - Xi(0) ] 

i=m \0<t<S-*K J 

= E P (f^m > ^Klog (j^r—) - s-'k 



K-i + l 

i=m 

i=m \ \ / 4 



i=l 

Note that we could drop the positive part in the last identity due to our assumption m = 2K/3 + 1. 
Now, we use the fact that for i < K - (2K/3 + 1) + 1, one has \og(K/i) - 1 > log 3 - 1 > .09, 
together with the second inequality in (III. 49) to obtain 

P (&K/S+1 < S- 2 K) < ^ * (-09^ < ^ e -0081K/4 < llv ^ e -*/500 

This finishes the proof of the Lemma. □ 

To complete the proof of Theorem we recall from |37j that a weak solution for the Atlas 
model as in Theorem |1.1| can be obtained by the following application of Girsanov's Theorem. 
Let Z±, . . . , Zk be independent Brownian motions such that Zj(0) = Xj(0), i = 1,...,K. Set 
T = 5~ 2 K and let Q° denote their joint law during the time interval [0, T] on the canonical sample 
space of continuous M^-valued functions on [0, T] with the usual Brownian filtration. Consider the 
martingale 

K ,t 



M(t) =J2f hZi{s)=z w {s)}dZi{s), t > 
»=i Jo 
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Note that its quadratic variation at any fixed time t > is given by 

(M) (t) 



(j2 f Q 1 {z 1 (s)=z w (s)} ) j ds = t, 



since two independent Brownian particles can simultaneously be the leftmost ones only on a set 
of Lebesgue measure zero. Hence, by Levy's characterization of Brownian motion, the process M 
is a standard Brownian motion under Q°. If we now change the measure to Q s according to the 
formula 

^ = exp (SM(T) - ^ , (111-50) 

then, under Q s , the law of the process (Z\(t), . . . ,Zx(t)), t G [0, T] is that of the Atlas model 
during the time interval [0, T]. 

Let J~t denote the cr-algebra generated by the entire process (Zi(t), . . . , Zx(t)), t G [0, T\. 
Clearly, one has the decomposition 

J~T = Qt V T~Lt, 

where Qt is the cr-algebra generated by the top 5 J indexed coordinate processes 

(Z K - 5J+1 (t),...,Z K (t)), t G [0, T], 

%t is the cr-algebra generated by the rest of the coordinate processes, and V refers to the smallest 
cr-algebra containing the two. 

By our assumption, 5J < K/3. Hence, the process (Zi(t), . . . , ^2Jr/3+l(i))> t £ [0>^1 i s mea- 
surable with respect to %t- For any fixed t G [0,T] define 

Z(-[\(t) = min Zj(t). 
(1JV ' i=l,...,2A73+l 

Then, the process Zr\\(t), t G [0,T] is also measurable with respect to Ht- 

Now, consider an arbitrary ^-measurable function F such that < F < 1. To simplify the 
notation, we will denote expectations with respect to the measures Q and Q s by Q(-) and Q S (-), 
respectively. By Lemma 3.3 and the change of measure formula (III. 50) we have 

Q s ( F ) < Q s ( f M* w >t } ) + s (i{* 2 k /3+1 <t }/ 

. . (111. ol J 

< go (Fe SM ^- s2T /\, 2K/w>T} ) + llVKe-*/™. 

Now, on the set {&2K/3+1 > the process M(t), t G [0,T] is identical to the process M(t), 
t G [0,T], where the latter is defined by 

2K/3+1 t 

M(t)= Yl / \ Zi (s)=z mi s)} dZ i(^ te[0,T}. 
i=i - 70 

Hence, it holds 

Q° (^ M(T) - 52r/2 l { . 2 K /3+1 >T } ) = Q° (^^-^/ 2 1 { , 2K/3+1 >t } ) < Q° (Ve^)"^ 

Note that A1(T) is measurable with respect to T~Lt while F is measurable with respect to 
Qt- Moreover, under Q°, the cr-algebras Qt and Ht are independent of each other. Using this 
observation and \M \ (T) = T, we obtain 



Q0 ^ Fe SA{T)-PT/2\ = Q° {F) Q° ^(^) = Q0 (i?)> 
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Combining this with the previous inequality and (III. 51), we get 



Q 5 (F) < Q°(F) + llVKe~ K/50 °. 



(111.52) 



For the rest of the argument we will assume that the particles indexed by the top 5 J indices are 
independent Brownian motions starting from their respective initial conditions, the idea being that 
all probabilities under the actual measure Q s can be bounded from above as in (III. 52). 



Back to the K particles Atlas model, consider the event {04J+1 > 5~ 2 K} as in Lemma 3.2 
On this event, during the time interval [0,(5 -2 -fr] the top J processes are identical to the top J 
processes among those that started at the top 5 J positions at time zero. Let Y(, . . . ,Y^j be the 
ranked processes Xk—5J+i, ■ ■ ■ ,Xk hi the increasing order. Also, set 

a (t) = , t > 0. 

Ei=2 lo g * 

Then, by Lemma |3.2| 

P (a(t) = a'(t), for all < t < 5~ 2 K) > 1 - 6 J 2 K~ 3/2 e- K/5 ° . (111.53) 

We now prove a concentration of measure property for a'(t), t G [0, 5~ 2 K\. Relying on ( III.52[ ), 
we can assume first that Xx-5j+i, ■ ■ ■ ,Xk evolve according to independent standard Brownian 
motions. If we let 

(£(*): i = l,...,J-l) = (Y>j_ t+1 (t)-Y^ J _ l (t): i = 1, . . . , J - 1), te%8~ 2 K}, 

then (£'i(t), . . . ,£j_i(t)), t £ [0, 5~ 2 K] can be viewed as a vector of (J — 1) component processes 
of a reflected Brownian motion in the (5 J)-dimensional positive orthant with zero drift vector and 
a constant diffusion matrix. 

By (III. 45) and the paragraph following it, we know that the process £'(£), t G [0, 5~ 2 K] 
satisfies a QTCI with respect to the norm 11-11,5-2^ max with the constant fj,'J s 5~ 2 K. Here, we have 
abbreviated 5 3 x 400 by //. 

Note that for any fixed t G [0, 5~ 2 K], a'(t) can be written in terms of as 



a'(t) 



Hence, the function that takes the paths of £'(t), t G [0,5 2 K] to the paths of ct'(t), t G [0,5 2 K] 
is Lipschitz with respect to ||- 1^-2^- max norms with the Lipschitz constant 



log( J!/*l 



C a (J) 



It follows that the random variable a! := sup 0<t< g-2^ a'(t) can be viewed as the image of £'(t), 
t G [0, 5~ 2 K] under a Lipschitz function with the Lipschitz constant C a (J). Thus, by Lemma 
its law satisfies a QTCI with the constant 



2.2 



C a (J) 2 n'J 3 5- 2 K. 

Let m a be the median of a' under Q°. Setting \i = 8/x', we deduce from Theorem 



2.1 



and 



(111.52): 



Q 5 (a' > m a + rV^j < exp ^~ 



r 2 5 2 
V<Sa(J) 2 J 3 



+ llVKe- K / 500 
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for all r greater than a constant depending only on J and 5. Combining this estimate with (III. 53) 
(and bounding J by K) we get 

Q 5 (a>m a + rVK) < exp ( ) + ll^e"^ 500 + Q^Ke~ K ^. 

v 1 V fJ,C a (J) 2 J 3 J 

The observation that the sum of the last two summands is smaller than the first summand for all 
sufficiently large K yields 

/ \ / r 2 5 2 \ 
Q s [a > m a + rVK) <2exp = . 



We note that m a is not the median of a under Q s . However, by (III. 52) and (III. 53) one has 

Q 5 (a < m a ) <Q°{a< m a ) + 11 VKe"^ 500 

< Q° (3' < m a ) + llVKe- K / 500 + Q^Ke~ K ^ 

< l/2 + l7^/Ke- K/500 . 



This completes the proof of Theorem |1.1| □ 

4 The infinite rank-based system 



This section is devoted to the proof of Theorem |1.3| The first step in the proof is to understand the 
dynamics of the process (Xn\(t), . . . ,X( n )(t)), t G [0, T] of the n leftmost particles in the particle 
system of Theorem |1.3| 



Lemma 4.1 There exist stopping times = tn < t jv+i < • • • such that the following is true. 

(a) linim^oo t N+m = oo with probability one. 

(b) For each m G N there exists a system of i.i.d. standard Brownian motions /3^ m \ . . . , such 
that for all i = 1, . . . ,n one has the dynamics 

dX(i)(t A T N+m ) = l{ Tjv+m >t}^ dt + dp\ m \t A T N+m ) + -dL^ x ,i){t A T N+m ) 

-^dL{i,a-i){t A T N+m ) 

on [0, T\. Hereby, a A 6 denotes min(a, b) for any two real numbers a, b. 

Proof. We define inductively the sets Ajv C Ajy +1 c . . . and the stopping times = r^v < tjv + i < 
• • • by 

AiV+ m = {k>l\31<i<N, 0<s<T N+m : X k (s) = X {l) (s)}, (IV.54) 
Tw+m+i = inf{s > T N+m \ 3 1 < i < N, k £ A N+m : X k (s) = X^(s)} (IV.55) 

for all m = 0, 1, 

The proof of Proposition 3.1 in [42j shows that, with probability one, it holds lirrim^oo TN+ m = 
oo and the sets A]y +m are finite for all m = 0, 1, . . . . Moreover, the same proof implies that for each 
such number m the paths of the process X^(t A rjv +m ), . . . ,Xr n )(t A TN +m ), t G [0, T] are given 



by the paths of the n leftmost particles in a particle system as in (1.1) with / = {!,..., |Ajv +m |} 
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which is stopped at time T^ +m . Hence, by Lemma 4 in |37j we conclude that assertion (b) of the 
lemma is true for our choice of the stopping times = t/v < tjv+i < □ 



Next, fix a K G N and let XLy . . . ,X'^ K ^ be the ranked particles in the system (1.1) with 



I = {!,..., K}. Also, let L' 



L' 



(1,2)' • • -^(K-1,K) 



be the local time processes at zero of the spacings 



processes in that system. Recall the definition of the norm ||-||t,2 i n ( 1-9 ) . From Lemma 
deduce the following concentration of measure property of the finite particle system. 



2.4 



we can 



Corollary 4.2 Let A and B be measurable subsets of (C([0,T],M n ), ||.||r,2) such that 



[ 1>2) (t),...,L[ n _ ltn) (t)), te[Q,T\)€A) >-, 
p(((l[ 2) (t)-l| 1) (t),...,l ( ' n) (t)-l( n _ 1) (t)),ie[0 1 T])eB) > l - 

and for any r > set 



A r = {he C([0,T],I 
B r = {he C([0,T};. 



on— In 



on— !■* 



jnf \\h - W\t,2 < r}, 
jnf \\h — W\t,2 < r }- 

heB 



(IV.56) 
(IV.57) 

(IV.58) 
(IV.59) 



Then for all n > 2 4 ^^gp^ 2 and r 2 > 2 4 3^ 2( ^T'° g2 « holds 
p(((L' (li2) (t), . . . , L[ n _ 1>n) (t)), t G [0, T]) ^ An) < exp ( - ^~^ T 
P(((x ( ' 2) (t) - . . . ,X{ n) (t) - X'^t)), t G [0,21) i B r2 ) < exp ( - ~ ^ 

Proof. From the considerations in step 1 of the proof of Theorem |1.1| we see that the process 

(L{ li2) (t),...,L' (n _ 1)n) (t)), te[0,T] (IV.60) 

can be obtained by applying the map <3?£ ( ' aR d then the canonical projection of C([0, T], R^ -1 ) 
onto C([0, T], W 1 - 1 ) to the process 

(5 2 -S x )t + (3' 2 (t) -^(t),...,(5 K - 5 K ^)t + (3' K (t) - , t € [0, T] 



(IV.61) 

Hereby, the Brownian motions /?£,..., are defined analogously to the Brownian motions f}\ , . . . , /3k 
in the proof of Theorem There, we have seen that the process in (IV.61 ) satisfies a QTCI with 
respect to the norm ||.||t2 with the constant 16(K — 1) _1 T. Combining Lemma 2.2 and Lemma 



2.4 we conclude that a QTCI with the constant 



K - 1\2 

n 



1)1 



2 16T _ 2 6 (K - ifT 
' K-l ~ n- 1 



(IV.62) 



2.1 



we 



applies to the process in (IV.60) with respect to the norm ||.||t,2- Hence, from Theorem 
obtain the first of the two claimed concentration of measure results. 

To show the second one, we recall from step 1 of the proof of Theorem |1 . 1 1 that the process 

(X{ 2) {t) - X' {1) (t), . . .,X[ n) (t) - X ( ' n _i)(i)), t G [0,T] (IV.63) 
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is the image of the process in (IV. 61 ) under the successive application of the map and the 

— W^ 1 ) onto CQO.T^R 71 - 1 ). 



canonical projection of C([0,T 



Moreover, we can rewrite (11.17) as 

(^ C "" 1) (/ l ))(t) = h(t) + (/<*-!) - Q^- 1 ))(($f X_1) (/»))(*)), t G [0,21 (IV.64) 

for all ft G C([0,T],M^- 1 ), where I^" 1 ) is the (K — 1) X (K — 1) identity matrix. It follows that 
for all fti,ft 2 G C([0, T], one has the estimates 

)(K-1) 

//i )) 

T,2 



(*g ( * _1) (^)) 



< ||ft 2 - fti||r,2 + 



< ||ft2 - ftl||T,2 + 



(J<*-1) _Q(^-1))(($^ - -(ft 2 ))(. 



3V2 



(ft 2 )) - ($ 



(fcl)) 



T,2 



(M)(0) 



T,2 



< f 1 + 3\/2(ii: - 1)1) ||ft 2 - fta|| Tj2 . 



In the second inequality we have combined the fact that the matrix I^ K ~^ — Q( K ~^ is tridiagonal 
with the elementary inequality («i +a 2 + 03) 2 < 3(a 2 + a 2 01, «2 ; a 3 G K. The third inequality 

is a consequence of Lemma |2.4| 



Now, it follows from Lemma 2.2 that the process in (IV. 63) satisfies a QTCI with respect to 



the norm ||.||t2 with the constant 

2 

\ + 3V2(K -1) 



K-l 



n — 1 



162 < 2 6 3 2 (K - 1) 5 T 



K-l 



n 



1 



The second claim of the corollary is a consequence of this and Theorem 2.1 



(IV.65) 
□ 



The last ingredient in the proof of Theorem 1.3 is an estimate on how fast the stopping times 
= T/v < Tjv+i < • • • in the proof of Lemma 4.1 grow to infinity in terms of the initial positions of 
the particles. 

Lemma 4.3 Let the Assumption 1.1 be satisfied with a constant c > 0. Then for all natural 



numbers m > 



maxj=i,...,M-i \Sj-Sm\T 



W(T N+m <T)< 



+ 1 (= 
NT 



AT 



+ 1) one has the inequality 



c(cm — c — AT) 



■ exp 



1 
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cm 



AT)' 



(IV.66) 



In particular, there exists a constant C(c, M, n, T, A) > independent of m such that 



<T)< C(c,M,n,T,A)e- 



m G N. 



(IV.67) 



Proof. We fix a natural number m as in the first statement of the lemma and note that on the 
event {r^v+m < 2} there exist numbers 1 < i < N and j > N + m such that the particle, which 
was the i-th from the left in the initial particle configuration, appears on the right or at the same 
position as the particle, which was the j-th from the left in the initial particle configuration, at a 
time t G [0, T]. Using this observation, the union bound and the definition of A (see the statement 
of Theorem 1.3), one has the chain of inequalites 

N 00 

Hr N+m < T) < p ( su p (^ft) " w &)) >~ AT + X M ~ X M) 



i=l j=N+m 



0<t<T 



< N Y F( sup (Wi(i) - Wj(t)) > -AT + X,(0) - X N (0)). 



j=N+m 



0<t<T 
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From Bernstein's inequality for Brownian motion (see page 145 in |38j). Assumption 1.2 and 
Ac 



the assumption m > ijr + 1 it follows that the latter expression can be bounded further by 



oo 

N ex p( 



j=N+m 
oo 



(*,(())- XaKQ) -Ar)- 

2T 



< 



oo 

J] exp ( 



k=m 



(cfc - Ar) s 

2T 



iV^exp(-— (A; 



Next, we note that the latter integral is equal to the probability that a standard normal random 

m-i-^ / — r 
variable exceeds ^= c multiplied by W27r^. Using this and the standard estimate 



f°° _£i , 1 _Z 

/ e 2 dz < -e 2 , y > 



(IV.68) 



one ends up with the first statement of the lemma. 



Finally, to see (IV.67), it suffices to observe that the argument of the exponential function on 



the right-hand side of the inequality (IV. 66) is a quadratic polynomial in m, in which the coefficient 
of m 2 is given by □ 



We can now prove the following refined version of Theorem 1.3 



Proposition 4.4 Let the sets A, B, A r , r > and B r , r > be defined as in Theorem 1.3 and 



let the constant C(c, M,n,T, A) be as in Lemma 4-3, Moreover, let m\,m2 € N be such that for 

_ Z 2 

all natural numbers m > mi (or m > 7x12) the value of C(c, M,n,T, A)e 3T m is less or equal 

to the difference between the left-hand and the right-hand side of the inequality ( I.lCty (or (1. 11), 

respectively). Also, define C\ and C2 to be the smallest positive real numbers such that 



3(n - 1) 
2 9 c 2 



C{ > mi, 



Vr > Ci : r > 2 Z 



4 _ /2(N + mi(r) - l) 5 Tlog2 
n — 1 



n 



l\r 2 

C 2 7 > 7712, 



2 9 3c 2 , 
r > C 2 : r > 2 4 



2(JV + mi(r) - l) 5 Tlog2 
n — 1 



(IV.69) 
(IV.70) 
(IV.71) 
(IV.72) 



where we have set mi(r) = ( ^ggr^ J ? r7 anc ^ m<2( r ) = ( 993c 2 ) 7 r7 • i/iere exisi positive 

constants C3, C4 depending on c, A, M, n, T and the value on the left-hand side of (1. 10) and 



(1. 11), respectively, such that for all n > C\ and r2 > C2 it holds 

((£(1,2) (t), • • • , L(„-i,„)(t)), t 6 [0, T]) £ 4 n ) < C 3 exp 
((X (2) (t)-X (1) (t),...,X (n) (t)-X (n _ 1) (t)),ie [0,T]) £B r2 ) <C 4 exp 



.2 10 
I (n — 1)7 c 7 



!8 5 „ 

2^3rT 



2 10 



1 (n — 1) 7c 7 

2 ' 18 9 ~~ 

2T37T 
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Proof. Since the way of proof is the same for both inequalities, we only provide the proof of the 
first one. To this end, we fix an n as in the statement of the proposition and set 



mi 



3{n-l) 
29c 2 
N + mi. 



i > 



(IV.73) 
(IV.74) 



We assume from now on that r\ is such that rh\ is an integer. If this is not the case, one merely 
needs to replace m\ by the smallest integer which is larger than rh\. Moreover, we note that the 
inequalities r\ > C\ and (IV. 69) imply rh\ > m\. 

Next, we recall the definition of the stopping time tk x an d observe 

{1>2) {t),...,L {n _ 1>n) (t)),te[0,T))£A r ^ 
<P(^((L ( i i2) (t),...,L (n _i jn) (t)),t G [0,T]) £A ri ,T Kl >TJ +F(TN+fh 1 <T) 
<p(((L ( i j2) (t),...,L ( „_i >n) (t)),t G [0,T]) ^A ri ,r Kl >t) +C(c,M,n,T,A)e-£™ 2 i, 



where the last inequality is a consequence of Lemma 4.3 



To obtain an upper bound on the first summand, which we call term (*), we remark first that 
for all triples of indices 1 < i < j < k it holds 



F(3t G [0,71 : X (t) (t) = X {j) (t) = X {k) (t)) = 0. 



(IV.75) 



Indeed, arguing as in the proof of Lemma |4.1| but replacing N by max(fc, M), we deduce the 
existence of stopping times = T max ( fc < 'r m ax(fc,Af)+i — • • • tending to infinity almost surely 
and such that for each I G N the dynamics of the k leftmost particles, stopped at T max (fc jv/) + /, is 
given by the dynamics of the k leftmost particles in a finite particle system as in ( |I.1[ ), stopped at 
7max(fc,M)+z- Hence, (IV.75 ) is a consequence of the considerations in section 2.2 of [26], Proposition 



1 in [26] and the fact that a countable union of P-null sets is P-null set. 

Next, we recall from the proof of Lemma |4.1| that for each m G N the paths of the process 
X(i)(t Arjv+m)) • • • ,X( n ^(t AT7v+m)> t G [0, T] are the paths of the n leftmost particles in a particle 
system as in (1.1) with | Ayv+ml particles, stopped at Tjv+ m - Moreover, (IV.75) shows that |A7v+m| 
Y ■ in Thus, on the event {tk x > T} the paths of the process 

(-£(1,2) (0i •••)£( n— l,n) 

(t)),tG[0,T] 

can be written as the composition of the map 1 ' with the canonical projection of C([0, T], M^ 1-1 ^ 
onto C([0,T],IR' 1 " 1 ) applied to the paths of the process 



{{52 - h)t + - Pr'it), (S Kl - S Kl -!)t + f3^>(t) - ^(t) j , t G [0, T], 

where Z?!" 11 ^, . . . , are i.i.d. standard Brownian motions defined in Lemma 

the proof of the first statement of Corollary 4.2 and using (IV. 70), we cone 



4.1 



Hence, following 



ude that term 



is 



bounded from above by exp 

All in 

above by 



rf(n-l) 



2 9 (Xi-l) 5 T / 

All in all. \\v have shown t hat P( ( (Xri )2 ) (*),•••, £( n -l,n) (*));* € [°> T D ^ A n) 



is bounded from 



exp 



r 2 (n — 1) 
2 9 (K 1 - IfT 



+ C(c,M,n,T,A)e-£ m *. 
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Plugging in the values of m\ and K\ one observes that the leading order term in the variable 

r\ is the same for both exponents on the right-hand side of the latter inequality and equals to 

4 a 10 

~ r \ ' is 5 — • This immediately finishes the proof of the proposition and completes the proof of 
Theorem 11,31 □ 
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